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Abstract Bats have increasingly been recognized as the nat- 
ural reservoir of severe acute respiratory syndrome (SARS), 
coronavirus, and other coronaviruses found in mammals. 
However, little research has been conducted on bat 
coronaviruses in South Korea. In this study, bat samples 
(332 oral swabs, 245 fecal samples, 38 urine samples, and 
57 bat carcasses) were collected at 33 natural bat habitat sites 
in South Korea. RT-PCR and sequencing were performed for 
specific coronavirus genes to identify the bat coronaviruses in 
different bat samples. Coronaviruses were detected in 2.7% 
(18/672) of the samples: 13 oral swabs from one species of the 
family Rhinolophidae, and four fecal samples and one carcass 
(intestine) from three species of the family Vespertiliodae. To 
determine the genetic relationships of the 18 sequences ob- 
tained in this study and previously known coronaviruses, the 
nucleotide sequences of a 392-nt region of the RNA- 
dependent RNA polymerase (RdRp) gene were analyzed phy- 
logenetically. Thirteen sequences belonging to SARS-like 
betacoronaviruses showed the highest nucleotide identity 
(97.1-99.7%) with Bat-CoV-JTMC15 reported in China. 
The other five sequences were most similar to MERS-like 
betacoronaviruses. Four nucleotide sequences displayed the 
highest identity (94.1-95.1%) with Bat-CoV-HKUS from 
Hong Kong. The one sequence from a carcass showed the 
highest nucleotide identity (99%) with Bat-CoV-SC2013 from 
China. These results suggest that careful surveillance of 
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coronaviruses from bats should be continued, because animal 
and human infections may result from the genetic variants 
present in bat coronavirus reservoirs. 


Keywords Coronavirus - Bats - South Korea - Severe acute 
respiratory syndrome - Middle East respiratory syndrome - 
Phylogenetic analysis 


Introduction 


Coronaviruses are large, pleomorphic, enveloped viruses con- 
taining a single linear, positive-sense single-stranded RNA 
molecule. Coronavirus genomes are approximately 27— 
32 kb in length, the largest continuous RNA genomes among 
mammalian viruses [1]. Coronaviruses are classified into four 
genera: alpha, beta, gamma, and delta [2]. Coronaviruses are 
the second most prevalent cause of the common cold in 
humans [3]. In livestock and poultry, coronaviruses are recog- 
nized causes of enteric and respiratory infections that are often 
fatal in young animals [4]. An outbreak of severe acute respi- 
ratory syndrome (SARS) in 2002 and 2003 resulted in infec- 
tion of 8096 people worldwide, with 774 (9.5%) of them 
dying from this novel human coronavirus [5]. SARS corona- 
virus is believed to have been acquired from an animal spe- 
cies, and the Chinese horseshoe bat may have been the orig- 
inal source of SARS infection [6-9]. The emergent Middle 
East respiratory syndrome coronavirus (MERS-CoV) may al- 
so have originated in bats, with fragments of MERS-CoV 
genes being identified in bats from both Saudi Arabia [10] 
and Africa [11]. More than 200 viruses have been isolated 
from and detected in bats [12]. The ability of bats to fly and 
migrate, as well as the large sizes of their social groups, pre- 
disposes bats to acquire and maintain viruses [13]. 
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Given the import of MERS into South Korea [14] and the 
presence of SARS in the relatively close geographic location 
of China [9] (Fig. 3), together with the fact that bats are a 
reservoir for coronaviruses, the prevalence of coronavirus in- 
fection in Korean bat species should provide valuable infor- 
mation. Although SARS-like and MERS-like bat 
coronaviruses have already been reported in the feces of 
Korean bats [14], the prevalence and presence of 
coronaviruses in other bat samples have not been described. 
Therefore, we investigated the prevalence of coronaviruses in 
Korean bat species using various kinds of samples, including 
oral swabs, fecal and urine samples, and carcasses. The 
coronaviruses found were characterized for their relationships 
to each other and to coronaviruses isolated in South Korea and 
in other countries. 


Materials and Methods 
Sample Collection 


A total of 672 samples (332 oral swabs, 245 fecal samples, 
38 urine samples, and 57 bat carcasses) were collected at 33 
sites of natural bat habitat from January 2016 to September 
2016 (Table 1). Bats were captured by net for collection of 
oral swabs, fecal samples, and urine samples, and released 
immediately after sampling. Fecal samples were also col- 
lected from guano. The bat carcasses collected were natu- 
rally dead bats found at the site where the other samples 
were taken. Swab samples were kept in viral transport me- 
dium at 4 °C, and organs collected from carcasses and other 
samples were stored at —80 °C before processing. 


Genome Detection and Sequencing 


All fecal and tissue samples were resuspended in 1% 
antibiotic-antimycotic solution (Corning, USA) diluted in 
phosphate-buffered saline (PBS), and clarified by centrifu- 
gation at 3500xg for 10 min. RNA from 200 uL sample was 
extracted with the QI[Aamp® Viral RNA Mini Kit (Qiagen, 
Germany) and eluted in 60 wL RNase-free water. cDNA was 
synthesized by PrimeScript First Strand cDNA Synthesis 
Kit (TaKaRa, Japan) following the manufacturer’s instruc- 
tions. Bat-CoV screening was performed by a pan- 
coronavirus PCR method based on primers used by Poon 
et al. (Corona 1 forward, 5'-GGTTGGGACTATCC 
TAAGTGTGA-3', and Corona 2 reverse, 5'-CCAT 
CATCAGATAGAATCATCATA-3’) [15], followed by se- 
quencing of the amplified product to confirm Bat-CoV iden- 
tification. The pan-coronavirus primers were used to ampli- 
fy and sequence a 440-bp segment of the highly conserved 
RNA-dependent RNA polymerase (RdRp) gene. Cycle se- 
quencing reactions were performed using the BigDye® 
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Terminator Cycle Sequencing kit version 1.1 (Applied 
Biosystems, Foster City, CA, USA). Reactions were puri- 
fied using SigmaSpin™ Post-reaction Clean-Up Columns 
(Sigma, St. Louis, MO, USA) and sequenced on an ABI 
PRISM 3130 automated capillary DNA sequencer 
(Applied Biosystems) according to the manufacturer’s in- 
structions. Primer sequences were excluded from the raw 
sequence data, and the sequences we analyzed were 392- 
bp long. All sequences obtained in this study were submit- 
ted to GenBank (accession numbers K Y432454—432471). 


Phylogenetic Analysis 


The nucleotide sequences were aligned and compared to 
52 selected human and animal CoV sequences available 
from the GenBank database using ClustalW software im- 
plemented in BioEdit version 7.0.9.0. The phylogenetic 
trees were drawn using the neighbor joining method using 
the maximum composite likelihood model with MEGA 7 
software. 


Results 
Prevalence of Coronavirus Infection in Bats 


Coronaviruses were detected in 13 oral swabs, one carcass 
(intestine), and four fecal samples from a total of 672 sam- 
ples, giving an overall detection rate of 2.7% (18/672). 
Four of the 14 bat species tested were found to harbor 
coronavirus. Thirteen positive oral swabs were detected 
from one species of the family Rhinolophidae, which 
accounted for 24.3% (163/672) of the samples, and the 
other positive samples were from Vespertilionidae bats that 
comprised 55.5% (373/672) of the samples. 

The 13 coronavirus-positive oral swabs were collected 
from bats at an abandoned mine in Jeonbuk province (Fig. 
1). Oral swabs and other samples (7 = 60) were obtained from 
three species of bats, Rhinolophus ferrumequinum, 
Miniopterus schreibersii, and Myotis macrodactylus, but 
coronaviruses were only detected in samples from 
R. ferrumequinum (Table 2). One carcass was collected from 
a cave in Chungbuk province (Fig. 1). In total, 94 samples 
were obtained from six bat species (M. schreibersii, Myotis 
petax, Vespertilio sinensis, Murina leucogaster, 
R. ferrumequinum, and M. macrodactylus), and coronavirus 
was detected from a V. sinensis sample. Four coronavirus- 
positive fecal samples were collected in four different regions, 
in Gyeongbuk province (Andong, Yeongdeok, and Gyeongju) 
and in the metropolitan city Gwangju (Fig. 1). In the forest of 
Andong, 43 samples in total were collected from three species 
of bats (Eptesicus serotinus, M. petax, and Pipistrelus 
abramus), and coronavirus was detected from one 
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Table 1 Bat sample information 
Region | Region 2 Site name Site type Major bat species No. of samples 
Oralswab Feces Carcass Urine Total 
Gangwon Yeongwol YW | Forest Eptesicus serotinus 11 18 29 
Gyeongnam Changwon CW 1 Abandoned mine Rhinolophus ferrumequinum 7 
Gyeongbuk Gyeongju GJ 1 Abandoned mine Rhinolophus ferrumequinum 70 6 12 90 
GJ2 Forest Myotis aurascens 1 1 2 
GJ 3 In area of private houses Pipistrelus abramus 10 10 
Forest Pipistrelus abramus 5 5 
Andong AD 1 Forest Eptesicus serotinus 25 18 43 
AD2 Forest Eptesicus serotinus 2 2 10 
Yeongdeok YD 1 Abandoned mine Myotis macrodactylus 3 5 
YD2 Forest Pipistrelus abramus 5 5 
Yeongju YJ 1 Abandoned mine Myotis nattereri 23 6 1 30 
YJ2 Forest, stream Pipistrelus abramus 5 5 
Yeongcheon YC 1 Abandoned mine Miniopterus schreibersti 12 12 
Cheongsong CS 1 Forest, valley Myotis petax 4 1 5 
Gwangju Shinchon GwJ 1 Forest Eptesicus serotinus 11 6 1 18 
Seokgok GwJ 2 Cave Rhinolophus ferrumequinum 5 5 
Daegu Dalseong DS 1 Forest, stream Pipistrelus abramus 5 5 
Ulsan Ulju US 1 Abandoned mine Rhinolophus ferrumequinum — 12 2 8 22 
US 2 In area of private houses Eptesicus serotinus 1 1 2 
Forest Eptesicus serotinus 31 31 62 
US 3 Forest Myotis aurascens 9 9 
Jeonnam Gurye GR 1 Forest Eptesicus serotinus 12 12 
Muan MA 1 Abandoned mine Myotis macrodactylus 7 1 8 
MA U Cave Myotis formosus chofukusei 5 we 
Jindo JD 1 Abandoned mine Rhinolophus ferrumequinum 10 1 1 12 
Hampyeong HPU Cave Myotis formosus chofukusei 3 3 
Jeonbuk Sunchang SC 1 Abandoned mine Rhinolophus ferrumequinum 43 3 10 4 60 
Chungnam — Seosan SS U Abandoned mine Myotis formosus chofukusei 3 4 7 
Chungbuk Danyang DY 1 Cave Miniopterus schreibersii 7 65 22 94 
DY 2 Cave Myotis macrodactylus 2 2 4 
DY 3 Cave Miniopterus schreibersii 18 1 6 25 
DY 4 Cave Miniopterus schreibersii 25 25 
DY 5 Forest, valley Myotis ikonnikovi 1 
DY 6 Forest, valley Myotis ikonnikovi 8 6 2 16 
DY7 Forest Murina leucogaster 7 10 17 
Total 332 245 57 38 672 


P. abramus bat. In Yeongdeok, a total of five samples was 
collected from P. abramus living in the forest, and one sample 
was positive for coronavirus. In the forest of Gyeongju, five 
samples in total were collected only from P. abramus bats, and 
one of the samples was positive for coronavirus. In the forest 
of Gwangju, 18 samples in total were collected from two 
species of bats (E. serotinus and Myotis aurascens) and coro- 
navirus was detected from one E. serotinus bat (Table 2). 


Phylogenetic Characterization of Bat Coronaviruses 
Isolated in Korea 


To determine the genetic relationships between the 18 bat 
coronaviruses obtained in this study and previously known 
coronaviruses, the 392-nt RdRp sequences were analyzed 
phylogenetically. All the sequences obtained in this study 
belonged to the betacoronavirus genus (Fig. 2). The sequences 


a Springer 


Lee S. et al. 


(%L'T) TLI/B1 
6/0 

(%V'ID 6/1 
(%L'L) 6€/€ 
£/0 

L7/0 

67/0 

vec/0 

vI/0 

8/0 

1/0 

€/0 

£c/0 

(%9°0) €9T/1 
(%0'8) €91/€1 
Lc1/0 


[B10], 


(%$°0) T8I/I 


A%UTL 6/1 


S1/0 


v/0 
v1/0 


€/0 
S/O 


0c/0 
cOT/0 


yngsunya 


L/0 
v/0 


€/0 


weusunys 


(%L'I7) O9/ET 


L1/0 


qT VS) veel 
61/0 


ynquoar 


cr/0 = S6/0 
S/0 
c/0 
v/0 
6/0 
S/0 
1/0 L/0 
cl/0—c9/0 
0/0 


wreuuosr = ues[Q, 


S/0 


s/0 


nsovq 


(%E'r) €7/1 


c/0 


1/0 


2(%6'S) LU/I 
£/0 


nfsuemy 


(%v'V) C/E 


ASUID LUE 
1/0 

8/0 

67/0 

v/0 


c/0 


81/0 
8v/0 
6L/0 

9/0 


L/0 


L/0 


67/0 


s/0 


vc/0 


ynqsuoskH ~=weusu0sky = uoMsuURD 


(QUISOyUT) ssvoIeD , 


uMmouyuy, 
sisuauls outiadsa, 


snupign snja1sidig 
SNIUDYISHD Snyassidig 
xpjad syody 

Liadayjou SOA 
SNIAJIDPOIDUL SOA, 
1A0yIUUOY] SOA 
lasnynfoys snsoutsof syodyy 
SUJISDAND SOA 
SISUALINSSN DULINA, 
daJSDBOINA] DULINAL 
snuyosas snaisajdy 
wunuimnbaunsiaf snydojouryy 


MSdaqia4yos snsajdoluyy 


satoods 


SQEMS [PIO g 


so]dures [e901 , 


[R10], 
uMOUUA) 


sepruoriodsa,, 
oeprydojouryy 
sepLoidoranyal 


Aqruuey 


sono ueyfodonouw/ssoutAold pure soioeds yeq JUSJOYIP Ul UOINGLYSIP snMARUOIOD = Z BI qu, 


pringer 


Qs 


Genetic Characteristics of Coronaviruses from Korean Bats in 2016 


™~ 
A Gangwon 
~ ‘ 
Seoul? 
Gyeonggi j j 
“= Chungbuk = 4 
Chungnam 
° rin 4 
Daejeon 
v~ Gyeongbuk 
- A 
D 
Jeonbuk é ee} A 
Ulsan 
ca ) Gyeongnam 
Gwangju ‘ : de 
A 
Jeonnam 


Jeju 


Fig. 1 Location of South Korea and geographical location of coronavirus-positive bat samples. Black circle Bat-CoV from oral swabsm. Black up- 
pointing triangle Bat-CoV from fecal samples. White up-pointing triangle Bat-CoV from a bat carcass 


from 13 oral swabs (16BO122, 129, 132, 133, 134, 135, 139, 
141, 143, 145, 150, 154, and 155) were grouped with the 
SARS-related bat betacoronaviruses, lineage B. Sequences 
from a carcass (16BT3) and four fecal samples (16BF104, 
109, 211, and 244) were grouped with MERS-like bat 
betacoronaviruses, lineage C (Fig. 2). 

The 13 sequences from oral swabs were 97.1—100.0% 
identical with each other, and the closest known CoV strain, 
with 97.1-99.7% nucleotide identity, was Bat-CoV-JTMC15 
that was recently isolated from R. ferrumequinum in China. 
The nucleotide identity with SARS-CoV was only 86.7— 


89.0%. One sequence from a carcass showed the highest iden- 
tity (99%) with Bat-CoV-SC2013 that was isolated from 
Vespertilio superans (synonym of V. sinensis), also from 
China [16]. The nucleotide identity with MERS-CoV was 
83.9%. Four sequences from fecal samples shared 96.9— 
100.0% identity with each other and showed 94.1—-95.1% 
identity with Bat-CoV HKUS5 [17]. The nucleotide identity 
with MERS-CoV was 84.1-84.4%. 

Due to shorter reported sequences from Kim et al. [14], 
another phylogenetic tree was constructed from the 355-bp 
RdRp sequences of Korean bat coronaviruses obtained by 
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Fig. 2 Phylogenetic tree 
constructed from 392-nt RdRp 
sequences of Korean bat 
coronaviruses detected in this 
study and from other 
coronaviruses. Phylogenetic trees 
were constructed by using the 
neighbor joining method and 
bootstrap values were determined 
by 3000 replicates. Scale bar the 
estimated genetic distance of 
these viruses. Black circle se- 
quences from oral swabs Black 
up-pointing triangle sequences 
from fecal samples. White up- 
pointing triangle sequences from 
carcass 
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Kim et al. [14] and in this study to analyze the genetic relation 
between Bat-CoVs in South Korea (Fig. 3). Thirteen se- 
quences from oral swabs were clustered with Bat-CoV B15— 
21, which was detected in fecal bat samples collected from an 
abandoned mine in Gangwon province. The major bat species 
present there was R. ferrumequinum, but the coronavirus- 
positive bat species were not specified [14]. Nucleotide iden- 
tity between Bat-CoV B15—21 and the 13 clustered sequences 
was 97.7—100.0%. The other three Korean bat coronaviruses 
(Bat-CoV B15-8, B15—40, and B15—41) belonged to the 
alphacoronavirus genus, but we did not detect coronaviruses 
belonging to this group in our study. 
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Discussion 


There are approximately 24 species of bats in South Korea, 
and most are insectivorous and relatively small [18]. Bats in 
South Korea usually do not migrate over distances greater 
than 1000 km [19], but apart from two species found only in 
South Korea, they are generally distributed in neighboring 
countries such as China, East Siberia, Japan, and Taiwan. 
Although it is not clear if bats go between China and other 
countries, some species that migrate over 500 km can carry 
pathogens to South Korea. The fact that the 18 bat 
coronaviruses (Bat-CoVs) detected in this study were all 
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Fig. 3. Phylogenetic tree of Korean bat coronaviruses detected by Kim 
et al. [13] and in this study, based on 355-nt sequences. Phylogenetic trees 
were constructed by using the neighbor joining method and bootstrap 
values were determined by 3000 replicates. Scale bar the estimated 


closely related to the Bat-CoVs identified from bats in China 
suggests the possibility that bat coronaviruses have been in- 
troduced from China and other countries into South Korea, or 
vice versa. This also can indicate the migration of bats over 
borders between countries. 

In this study, the prevalence of coronavirus in Korean bat 
samples was 2.7% overall. The prevalence by sample type was 
5.4% (13/332) for oral swabs, 1.8% (1/57) for carcasses, 1.6% 
(4/245) for fecal samples, and 0.0% (0/38) for urine samples. 
Reported rates from other countries for coronavirus prevalence 
in bats were 6.5% (64/985) and 5.3% (50/951) for China, 8.2% 
(16/195) for Italy, and 5.3% (32/606) for Mexico [20-23], not 
greatly different from our figure of 2.7% for prevalence in 
South Korea. 

Fourteen out of 24 Korean bat species were screened for 
CoVs in this study. The species positive for Bat-CoVs were 
R. ferrumequinum, V. sinensis, P- abramus, and E. serotinus. 
The Bat-CoVs that showed the highest nucleotide identity with 
our isolates were obtained from the same species that we iden- 
tified as positive, except for E. serotinus. To our knowledge, 
this is the first time a Bat-CoV has been detected from 
E. serotinus in Asia. In Europe (Italy), a lineage C 
betacoronavirus was detected from £. serotinus in 2014 [24]. 
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Bat-CoVs were usually detected at higher rates in alimentary 
samples (fecal samples, rectal swabs, and intestine) than in 
respiratory samples (oral and throat swabs) [7, 15, 17, 22, 25, 
26]. We collected alimentary and respiratory samples and urine 
samples. However, Bat-CoVs were not detected in urine sam- 
ples. Animal coronaviruses usually produce either enteric or 
respiratory infections. Animal models show similar clinical fea- 
tures to SARS and MERS patients; the majority of whom pre- 
sented with respiratory symptoms but in some cases also suf- 
fered from enteric complications [27]. Unlike other human co- 
ronavirus infections, a number of MERS cases were associated 
with acute renal failure. MERS-CoV replication in the human 
kidney suggests the potential of shedding viruses in urine [28]. 
Bat-CoV isolate Bat SL-CoV-WIV1 (KF367457) showed in- 
fectivity to Rhinolophus sinicus kidney cell lines, which means 
it may be possible to detect Bat-CoV in urine [6]. 

The current study revealed that Korean bats have SARS- 
and MERS-related Bat-CoVs, consistent with earlier work 
[14]. The previous study detected alpha- and 
betacoronaviruses only from fecal samples in four regions, 
two in Gangwon, one in Chungbuk, and one in Gyeongbuk 
[14]. The only betacoronavirus RdRp sequence Bat-CoV 
B15—21 (KU528590.1) was clustered with 13 Bat-CoVs 
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(16BO122, 129, 132, 133, 134, 135, 139, 141, 143, 145, 150, 
154, and 155) detected from Danyang, Chungbuk province. 
The identity between the 14 sequences was very high at 97.7— 
100%, even though the sample types and collection sites were 
different. The species of bat infected with Bat-CoV B15-21 
was not specified in the paper, but the major bat species of the 
sample collection site was R. ferrumequinum; the same spe- 
cies in which we detected 13 Bat-CoVs. The same Bat-CoVs 
were detected from fecal samples in the prior study [14] and 
from oral swabs in our study. Spike proteins, which are trans- 
lated from spike genes, define viral tropism by receptor spec- 
ificity and by membrane fusion activity during entry into cells 
[29]. Although the spike genes were not identified in this 
study, full spike gene analysis of Bat-CoV B15-21 suggested 
a low potential for human emergence. One MERS-like 
betacoronavirus, Bat-CoV B15-1-6 (KU528586.1, 
KU528585.1), was identified by Kim et al. [14], but the sim- 
ilarity of this isolate to our isolates could not be analyzed 
because the RdRp sequence was not obtained. High nucleo- 
tide identity between Bat-CoV B15-1-6 and Bat-CoV-SC2013 
(KJ473821) suggests that Bat-CoV B15-1-6 and 16BT3 might 
be very similar. 

Bat-CoVs have been identified in many countries including 
the Australasian region [10, 11, 14-16, 20-26, 30, 31]. But 
Bat-CoVs identified in Korea were genetically close to Bat- 
CoVs from China, a close neighbor of Korea. According to 
the results of Bat-CoV research in China, the degree of host 
restriction for coronaviruses in bat populations was high. For 
example, similar viruses were detected in the same species of 
bat in different regions, approximately 1600 km apart, while 
two different bat species in the same habitat had different 
coronaviruses [22]. In our results, three Bat-CoVs 
(16BF109, 211, and 244) from P abramus were clustered 
with each other, even though samples were collected from 
different habitats in different regions (Andong, Yeongdeok, 
and Gyeongju). The Bat-CoV detected from E. serotinus bat 
in the metropolitan city Gwangju (16BF104) was also clus- 
tered with these three Bat-CoVs mentioned above. However, 
there was no matching virus for 16BF104 found in the same 
species of bat in other locations, nor in other species in the 
region where samples were collected. Our data could support 
the possibility of relationships between genetic variation, geo- 
graphic locations where Bat-CoVs were isolated, and bat host 
species, but more data on Bat-CoVs in Korea would be need- 
ed to address this. 

In summary, the prevalence of coronavirus in Korean bats 
was determined and found to be comparable to or lower than 
the rates in other countries. We detected SARS-related and 
MERS-like Bat-CoVs that were genetically very similar to 
Bat-CoVs identified in China. For the study of Bat-CoV sea- 
sonal occurrence and transmission, continuous monitoring 
during all seasons is required, and ideally alimentary, respira- 
tory, and other samples need to be collected from each bat. 
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Future work to determine the complete genome sequences of 
Bat-CoVs from South Korea would give a more complete 
picture of their pathogenicity and the possibility of infectious 
spread to other animals and humans. 
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